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J> . Abstract. The complex WKB-Maslov method is used to consider an approach to the 

semiclassical integrability of the multidimensional Gross-Pitaevskii equation with an ex- 
ternal field and nonlocal nonlinearity previously developed by the authors. Although the 
WKB-Maslov method is approximate in essence, it leads to exact solution of the Gross- 
Pitaevskii equation with an external and a nonlocal quadratic potential. For this equation, 
If} ' an exact solution of the Cauchy problem is constructed in the class of trajectory concen- 

trated functions. A nonlinear evolution operator is found in explicit form and symmetry 
"*^h ' operators (mapping a solution of the equation into another solution) are obtained for the 

O j, equation under consideration. General constructions are illustrated by examples. 
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1 Introduction 



Experimental advances in the realization of Bose-Einstein condensation (BEC) in weakly inter- 
acting alkali-metal atomic gases have generated great interest in the theoretical study of the 
BEC. Its states and evolution are described using the Gross-Pitaevskii equation (GPE) [21IH] for 
the wave function fy(x,t) of the condensate confined by external field with potential V ex t(x, t) 
at zero temperature: 

-ihd t + ^ + V cxt {x,t) + x\^(x,t)\ 2 ^j $(x,t) = 0. (1) 

Here x G R™, p = -ihd/dx, t G R 1 , d t = d/dt, |^| 2 = 'F*^, 'F* is complex conjugate to <F, k is 
a real nonlinearity parameter, \^(x, t)\ 2 is the condensate density, and iV[\P] = J" Rri \ty(x, t)\ 2 dx 
is the number of condensate particles. Following quantum mechanics, we refer to the solutions 
of the GPE as states. 

Equation is of nonlinear Schrodinger equation type with the local cubic nonlinearity 
x\^/(x, t)\ 2 representing the boson interaction in the mean field approximation. Besides the 
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BEC, equation Q describes a wide spectrum of nonlinear phenomena such as instability of water 
waves, nonlinear modulation of collisionless plasma waves, optical pulse propagation in nonlinear 
media and others. In all these cases, space localized soliton-like solutions are of principal interest. 
However, the wave packets described by equation (J2) with V cxt = in multidimensional space 
(n > 1) with focusing nonlinearity (x < 0) are known to collapse, i.e. the situation where the 
wave amplitude increases extremely and becomes singular within a finite time or propagation 
distance (see, for example, Jl] for reviews). To eliminate the collapse, which is considered as an 
artifact of the theory, one has to consider some effects that would render collapse impossible. 
The nonlocal form of nonlinearity is significant since it can, basically, eliminate collapse in 
all physical dimensions (n = 2,3) [5]. Meanwhile, in the derivation of the GPE (^Q), a nonlocal 
nonlinearity term L n V(x, x! ) \^{x\t)\ 2 dx' V(x,t) arises, and it is reduced to x\^(x, t)\ 2 ^(x, t), 
x = fV(x)dx if the potential V(x,x') = V(x — x') is short-range [2J. 

In view of this consideration, of fundamental interest is the study of both the properties 
and the localized solutions of the nonlocal Gross-Pitaevskii equation. The integrability of the 
GPE is a nontrivial problem. In the one-dimensional case (n = 1) with V^xt = 0, equation ([T]) 
(called the nonlinear Schrodinger equation (NLSE)) is known to be exactly integrable by the 
Inverse Scattering Transform (1ST) method [HllZj- This is the only 1ST integrable case since 
with V ex t 7^ the 1ST fails even in the one-dimensional case. The same is true for both the 
nonlocal GPE in all dimensions and the local GPE Q in a multidimensional case (n > 1). 
Direct application of symmetry analysis [HI HE 1101 111! I12j to the nonlocal GPE is also hampered 
by presence of the nonlocal term and external potential in the equation. 

In |131 1141 1151 I16j a semiclassical integrability approach was developed for a generalized 
nonlocal GPE named there a Hartree type equation: 

{ -ihd t + H„(t)}y(x,t) = {-iWt + H(t) + xV(t,V(t))}y(x,t) = 0, (2) 

V(x,t) E L 2 (R%), V(t,*(t))= f dyy*(y,t)V(z,w,t)V(y,t). (3) 

Jr™ 

Here the linear operators ~H(t) = Tt(z, t) and V(z, w, t) are Weyl-ordered functions of time t 
and of noncommuting operators 

z = (p,x) = (—ihd/dx,x), w = (—ihd/dy,y), x, y £ M n , 



with commutators 



[z k ,Zj]_ = [w k ,Wj]_ = ihJ kj , [z k ,Wj]_ = 0, k,j = l,2n, (4) 

J = || J k j \\2nx2n is a identity symplectic matrix J = ( !jj 1 ) , I = I nxn is an identity 

(n x n)-matrix. This approach is based on the WKB-Maslov complex germ theory |18[ I19j and 
gives a formal solution of the Cauchy problem, asymptotic in formal small parameter h (h — > 0) 
accurate to 0(h N l 2S ) , where iV is any natural number. The Cauchy problem was considered in the 
Vt class of trajectory concentrated functions (TCFs) introduced in |13M14) . Being approximate 
one in essence, the semiclassical approach results in some cases in exact solutions. 

In the present work we construct an exact solution of the Cauchy problem in the Vt class for 
equation (j2j) with the linear operators 7i(z,t) and V(z,w,t) being quadratic in z, w: 

H(z,t) = ±(z,H zz (t)z) + (H z (t),z), (5) 
V(z,w,t) = -(z,W zz (t)z) + (z,W zw (t)w) + -(w,W ww (t)w). (6) 
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Here, TC zz (t), W zz (t), W zw (t), W ww (t) are 2n x 2n matrices, H z (t) is a 2n-vector; (•,•) is an 

n In 

Euclidean scalar product of vectors: (p,x) = Pj x j'-> P>£ G ( z , w ) = Yl z 3 w ji z i w e ^ 2n - 

By solving the Cauchy problem, we obtain a nonlinear evolution operator in explicit form. 
With the evolution operator obtained, we formulate a nonlinear superposition principle for the 
solutions of the nonlocal GPE in the class of TCFs. Also, we give symmetry operators in general 
form that map each solution of the GPE into another solution. The general constructions are 
illustrated by examples. 



2 Cauchy problem 

Here we give a brief account of the solution of the Cauchy problem for equations ©, ©, ©, 
and ©, following [HI HH EH1 

The class of trajectory concentrated functions 

To set the Cauchy problem for equations © , © , © , and © , following ^3 El EES] , we define 
a class of functions V\ via its generic element $>(x,t,ti): 



^(S(t,h) + (P(t,h),Ax)) 



(7) 



Here, the function <p(£,t,ti) belongs to the Schwartz space 5(R n ) with respect to £ G M. n , 
smoothly depends on t, and is regular in y/h for h — > 0, Af = x — X(t, h). The real function 
S(t,h) and the 2n-dimensional vector function Z(t,h) = (P(t,h),X(t,h)) define the V\ class, 
regularly depend on \fh in the neighborhood of h = 0, and are to be determined. The functions 
of the V\ class are normalizable with respect to the norm ||<£(i)|| 2 = (<fr(t)| <&(£)), where 

(tf(t)|*(*)) = f dxV*(x,t,h)$(x,t,h) (8) 

is a scalar product in the space At any time t £ I 1 the function $(x,t,h) G V l h is 

localized in the limit h — > in the neighborhood of a point of the phase curve z = Z(t, 0). For 
this reason we call V\ the class of trajectory concentrated functions. 
class transforms to the V\ 



For t = the Vt class transforms to the V2 class of functions 



^(x, K) = exp 1 1 [5(0, ft) + (Po(^), (a? - ^o(ft)))]} ^ ftj , 

(9) 



where Zq{K) = (Po(h), Xq(H)) is a point of the phase space R~™, and the constant S'o(ft) can be 
omitted without loss of generality. The Cauchy problem is formulated for equations ©, ©, 
©, and © in the V l n class of trajectory concentrated functions as 

^(x,t,h)\ t =o = ^(x,h), tp g v\. (10) 

Hamilton— Ehrenfest system 

To solve the Cauchy problem ©, ©, ©, ©, and we first obtain the Hamilton-Ehrenfest 
system (HES) of equations in the moments of the solution *(x, t, K) G V\ of equation ©. 
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The operators H(z, t), equation ©, and V(z, w, t), equation (jEJ), are self-adjoint, respectively, 
to the scalar product (JSJ) and to scalar product 



(*(t)|*(t)) R a» = / dxdyV*(x,y,t,h)<S>{x,y,t,h) 
in the space L 2 (M^). 

Define the mean value (A) for a linear operator A(t) = A(z,t) and a state ^?(x,t,h) as 
(A) = —^(*(t)\A\*(t)) = A*(t,h). 



From @, © we have 



where [A, £?]_ = AB—BA is the commutator of the operators ^4 and -B. We refer to equation (|12j) 
as the Ehrenfest equation for the operator A and function *$>(x, t, K) as in quantum mechanics [2D] • 

For A = 1, equation ((T2J) gives ||*(t)|| 2 = H^O)!) 2 = H^H 2 - This implies that the norm of 
a solution of equations ((21) and ® does not depend on time, and we can use the parameter 
k = x||^|| 2 instead of x in @. 

Assume that 

4 ,„ M) = »™, o6Zf , (13) 

are the moments of order \a\ of the function ^(x, t) centered with respect to z^(t, h) = (py(t, h), 
X\&(t,h)). Here {Az} a is an operator with a Weyl symbol (Az^) a , 

Azq, = z — zqt(t, K) = (Ap*, Ax*), Apq, = p — p^(t, h), Ax* = x — x*(i, h). 

Along with (|13j) we use the following notation for the variances of coordinates, momenta, and 
correlations: 
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A m (t,H) = \\A q , jk (t,h)\\ 2 nx2n = 7 rr i ^\\(^(t)\{Az j Az k + Az k Az j }\^(t))\\ 



211^"''' 

a pp (t,h) cr px (t,h) 
o- xp (t,fr) a xx (t,h) 

o- xp (t,h) = -\\({AxjAp k + Ap k Axj})\\ nXn , 

o- xx (t,h) = \\(AxjAx k )\\ nxn , o- pp (t,h) = || (ApjAp k )\\ nXn . (14) 

The Ehrenfest equations (|12j) in mean values can be obtained for the operators Zj, {Az} a and 
trajectory concentrated functions (J7| (see ^HE] for details). However, for solving the equations 
under consideration, equations ((2J), ©, ©, and ©, we need only equations for the first-order 
and second-order moments 

z 9 = J{H z (t) + [H zz {t) + HW zz (t) + W zw (t))]z*}, 

A* 2 = J[H zz (t) + xW zz (t)]A m - A m [U zz {t) + kW zz (t)]J. (15) 



Exact Solutions and Symmetry Operators 



5 



We call the system (|15|) the Hamilton-Ehrenfest system (HES) of the second order for equations 
J3J), (JHJ), and ©• Following ^2E], we equate the functional vector-parameter Z(t,h) of the 
V l n class of TCFs J7J) with zy(t,h), i.e. py(t,h) = P(t,h), x$,(t,h) = X(t,h). This relates the 
ansatz (J7J) to an exact solution of equation (J2J). 

Consider a phase space M, , dim A4 N = N = 3n + 2n 2 , of points g € .M with coordinates 

fl = (z,A)T, zgM 2 ™, z=(p,x) T , 

A = (A 4i ) T , Ay = Ajj, Ael n+2 " 2 , i,j = l,...,2n. 

Here A T is a matrix transposed to A. The coordinates of g G .M^ are written as matrix columns. 
The HES ()15|) can be considered a dynamic system in .M^: 

z = J{n z (t) + \n zz (t) + k(w zz (t) + w zw (t)M, (i6) 
A = J[n zz (t) + ^ 2Z (t)]A - A[n zz (t) + kw zz (t)]j. (17) 

With the substitution 

A(i) = A(t)A(0)A + (t), 
equation (|17|) is rewritten in equivalent form: 

i = J[H zz (t) + *W„(t)]A, A(0) = I (18) 

We call equation ()18|) a system in variations. 

Denote by fl(t,£) the general solution of equations (|T?t|) . (fTTj) : 

(t, £) = (ity, ft, £), ft, C), An (*, ft £), A 12 (t, ft, £), . . . , A 2n2n (t, ft, £)) T (19) 

and by g — the operator column 

g= (^,x,(Ap 1 ) 2 ,Ap l Ap 2 ,...,(Ax n ) 2 ) T - (20) 

Here 

£ = (Ci,...,Civ) T e M 3n+2n2 (21) 

are arbitrary constants. Given constants £, equation (JT5|l describes a trajectory of a point in 
the phase space M. N . 

Lemma 1. Let ty(x,t) be a partial solution of the GPE (J2J), (|3*|). (JSJ), © urei/i an initial solution 
^{%>~t)\ t — Q = ^(x). Define the constants (£(ty(t)) by the condition 

fl(t,C)= |j^pW)|fl|*(t)>, (22) 
and ifte constants £(?/>) by the condition 

(O,£) = -L(^>. (23) 



Taen, 

£(*(*)) = £(V), (24) 
i.e., £( v I / (t)) are integrals of motion for equations (J2J), ©, ©, and @. 
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Proof. By construction, the vector 

S(t) = |j^p<* (t)\mt)) = fl(*, £(*(*))) (25) 

is a partial solution of the system (|16[). (|17|). which coincides with g(t, £(?/>)) at the initial time 
t = 0. In view of uniqueness of the solution of the Cauchy problem for the system (|16(l and (|17[). 
the equality 

fl(t,etyO)=fl(*><W))). (26) 

is valid. The proof is complete. ■ 
Linear associated Schrodinger equation 

Let us substitute ©, © into (0), © and replace the operators z = (p,x) and w by A£ = 
5 - z(t, h, £) = (A# Ax) = (j?- P(i, ft, C), x - X(f, ft, <£)) and Aw = w - z(t, ft, £). Then we 
have 

{-t^t + «(t,*(t))}* = 0, (27) 
*(t)) = fl(*> *(*)) + (J5*(t, *(*))> M + \(&z^zz{t, *(t))Az), (28) 
fl(t, *(*)) = ^{Mt, H), [HzzW + *{W zz (t) + 2W w (i) + WWt))]z* (t, ft)) 

+ {U z {t),z^{t,h)) + -xSp(W^(i)A* 2 ), 

*(t)) = + (H zz (t) + xW„(t) + *rW„,,(t))z*(t, ft), (29) 
*(*))= W w (t) + &W„(t). (30) 

Let us replace the mean values zy(t,K), A^(i, ft) in the nonlinear GPE ((2*5|) by the respective 
terms of the general solution (|19j) of the system 1)16(1 . l)17j) . As a result, we obtain a linear 
equation: 

{-tftfit + £(t, s (t, £))}*(£, t, ft, £) = 0, (31) 
£(t, 0(t, £)) = 0(t, ft, <£) + ft C), Az) + i(Az, j5«(t)A5), (32) 

$j(t, ft, C) = ^(i, ft, £), [W„(t) + x(W^(t) + 2W^(i) + WW*))]^*, ^ £)) 

+ (W z (t), z(t, ft, £)) + ^xSp(W M (i)A(t, ^ £)), 

S*(t, ft, C) = W,(t) + (H zz {t) + *Wzz(*) + %W zw {t))z{t, ft, £), 

= H**(t) + ^W^(i), *(t, ft, £) = (P(t, ft, £), X(t, ft, £)). 

We call equation pT|) ifte linear associated Schrodinger equation (LASE) for the nonlinear Gross- 
Pitaevskii equation (|28|). More precisely, equation 1)31)1 is to be considered as a family of equa- 
tions parametrized by the constants £ of the form (|21() . Each element of the family 1)31)1 is 
a linear Schrodinger equation with a quadratic Hamiltonian with respect to operators of coor- 
dinates and momenta. Such an equation is well known to be solvable in explicit form (see, for 
example, |21) I22j). In particular, partial solutions can be found as Gaussian wave packets and 
a Fock basis of solutions, and Green function can be constructed. 



Exact Solutions and Symmetry Operators 



7 



LASE solutions and GPE solutions 

Consider a relationship between the solutions of the LASE and GPE. Let <&(x,t,h,<£) be the 
solution of the Cauchy problem for the LASE (|31|) 



t=0 



ip(x, h), 



(33) 



The function $(x,t,H,<£) depends on arbitrary parameters <£ which appear in the LASE (|31|) . 
Let <£ are subject to the condition (|23|) and are functionals £(V0- 



Theorem 1. The solution of the Cauchy problem Ijlfljl /or i/ie GP-E J2J is 



(34) 



Proof. The function <I>(x, t, K, £(ip)) satisfies equation (|31[) for arbitrary £ and also for (£ = £(V0- 
According to equation (|24|). in LemmaQwe do not violate the equality in (|31|) if we replace 
by £(*(*)). In view of equations (|23|) . (J2H1), one can see that $(f, h, £(?/>)) = $(x, i, K, £($>(t))) 
and the operator $)(t, g(t, <£)) in (|31|) becomes to $j(t, ^(t)) in (|27[). This implies that the function 
&(x,t,h,£(il))) satisfies equation (|27|) and the initial condition (|33|) . which correlates with 
Consequently, equation (|34j) is valid, and the theorem is proved. 

The relationship between the steps described above can be shown diagrammatically: 



{-ihdt + H„(t)}V(x,t) =0, 
*(x,0) = i>(x) 



Class T^of TSFs 



HES, general 
solution g(t,€) 




Family of LASE 

{-ihd t + h(t,£)}& = 0; 
Cauchy problem 
$(x,t,<t)\t=o =0(f) 



Theorem^makes it possible to obtain a nonlinear evolution operator for the GPE (|27|) in the 
7?£ class of TCFs ©. The evolut ion operator can be written as a nonlinear integral operator 
using the Green function of the LASE (|31|) with constants £ changed by according to 

relation l|3"4"Jl. 



3 Nonlinear evolution operator 

The Green function G^ix, y, t, s, g(t, <£), g(s, (£)) for the Cauchy problem (|31[). 1)33(1 is defined by 
the conditions 

H^+«(t,fl(t,C))]G a< (£,ir,t,s,fl(t,C),fl(a I C)) =0, (35) 
limG^(£,y,M,g(t,£),0(s,£)) =8{x-y). (36) 

Here the operator £j(t, £)), given by (|32|). is quadratic in coordinates and momenta. 
We shall seek for the required Green function under the simplifying assumption 

det£ pp (s) ^0. (37) 

Following, for example, |21| 122 j. we obtain 

GJx,y,t,s,8(t,£),g(s,£)) 



y/det{-i2irh\ 3 (t,s)) 
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xexp - 



S(t, h, Q (t, £)) - S(a, h, &(s, £)) + (P(t, h, <£), Ax) - (P(s, h, <£), Ay) 

1.(38) 



■i(Ay,Ai(t,s)A3 1 (t,s)Ay) + (Ax, Ag s)Ay) - -(Ax, A3 \t, s)X A (t, a) Ax) 



Here, Ay = y — X(s, H, £), nxn matrices A^(t, s), fc = 1,4, are blocks of the block matrix «4(£, s) 
of the system in variations: 

A=Jfi zz (t,h)A, A\ = = hnx2n, (39) 

^^(-1,1) <«> 

and 

s(t,H,B(t,£)) = J*{(P(t,h,£),x(t,h,e.)) -fj(t,nx)}dt. (41) 

Then, the following theorem is true: 

Theorem 2. Lei an operator U^it, s, •) aci on a given function ip(x), taken at an initial time s, 
as follows: 

U„(t, s, VO (x) = f G„(x, y, t, s, (t, Q(s, C(^))^(y) dV (42) 

ilere G H [x, y, t, s, g(t, C(V0, fl(s, £(^)))) is determined by (|38jl. (|41|). and i/ie parameters <L(ip) 
are obtained from 

B(t,£)\ t=s = 00(^ = ^(^0^). (43) 



Then, the function 

V(x,t) = U„(t,s,iP)(x) (44) 

is an exaci solution of the Cauchy problem for equations (|27|l. (|28[) u>ii/i i/ie initial condition 
^(x,t)\ t=s = ^(x), and the operator U^(t,s,-^ is the evolution operator for the GPE (|2"Tj) f/iai 
acis on i/ie cZass 0/ trajectory concentrated functions (J7J) . 

For the evolution operator 1)42(1 . the following properties can be verified by direct computation. 

Theorem 3. The operator U~ l (t, s, ■) , 

L>- 1 (t,s,^)(x)= ( G- 1 (x,y,t,s, (t,(^)),s(s,e:(v)))V(y)rf n y 

G„(x, y, s, t, S (s, £(^)), (t, (</>))) ^(y) d n y, (45) 

is i/ie Ze/£ inverse operator to U^(t, s, •) 0/ (|42|) . so t/iai 

i7- 1 (t,s,^(t, S ,^)(x)=^(x), tj) E V\. (46) 
Corollary 1. For a partial solution fy(x,t) of the GPE (|27|) . we have 

U„(t, s, (J" 1 (t, s, *(i))) (x) = *(x, t). (47) 
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Proof. Indeed, according to Theorem |2] we have ty(x,t) = U )t (t,s,'(p)(x), ip(x) = ^(x, t)\ t . 
In view of (|46jl . the left-hand side of equation (|47|) can be written as 

U„(t, s, U-\t, s, *(<)))(£) = U„[t, s, U-^t, s, U„(t, s, V>))] (x) = U„(t, s, iP)(x) = t). 

Thus, the statement is proven. ■ 

Theorem 4. The operators U^U,-) = U^CtjO,-) possess the group property 

U„(t + s,^)(x) = U„(t,y(s))(x), ¥(£,a) = Ck(s,^)(f). (48) 

Proof. The functions ty(x,t) = U^yt + s,tp^(x) and ty(x,t) of (|44[) are partial solutions of 
equation (|3T]) with the same trajectory Q(t,€) in the extended phase space, and is valid 
for the evolution operator of the linear equation (|31|) . Then, it is also valid for g(t, £(?/>)) 
corresponding to the nonlinear evolution operator (1421) , 
Substituting t + s — > t in (fl5|) , we find 

U x {t,il>)(&) = U x (t-8,*(8))(Z), (49) 
*(f,s) = ^(x) (50) 

s=0 

for a solution ^(x, t) of the Cauchy problem for the GPE (|27j) with the initial condition (|5()jl. ■ 



4 Symmetry and nonlinear superposition 

Symmetry operators 

A symmetry operator, by definition, maps a solution of an equation into another solution of this 
equation. Direct finding of symmetry operators for a given nonlinear equation is an intricate 
problem because of the nonlinearity of the determining equations. 

It is rare for this problem to be solved (see, for example, US])- The symmetry analysis 
of differential equations deals mainly with generators of one-parametric families of symmetry 
operators (symmetries of an equation) determined by linear equations [HI El 111! 1121 HUj . Using 
the evolution operator U^it, s,-) given by 1)42 Jl . we can formulate a general form for symmetry 
operators of the Gross-Pitaevskii equation (|27|). 

Let a be an operator acting in V\i (a : V\ — > V^) and ^(x,t) is an arbitrary function of the 
V\ class (%(x,t) G Vft). Consider an operator A(-), such that 

$(x,t) = A(*(i))(x) = U4t,W- 1 {t,^(t)))(x). (51) 

If ty(x,t) is a solution of the GPE (|27|) . then <&(x,t) is also a solution of equation (|27|). This 
follows immediately from Theorems El and [31 and Corollary ^ 

Thus, the operator A(-) determined by (|51|) is a symmetry operator for the GPE (|27|). 

Assume now that operator b and its operator exponent exp(ab) act in the V2 class, i.e., 
b : V\ -> P£ and exp(ab) : T 3 ^ ->■ Pg, where a is a parameter. 

Define a one-parametric family of operators B(a, •) via their action on an arbitrary function 

B(a, = Inexpiable/- 1 (^(i)))(x). (52) 

By analogy with the aforesaid, the operators B(a, •) constitute a one-parametric family of the 
symmetry operators of equation 1)27(1 . 
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It is easy to verify the group property 

B(a + /?, *(t))(f) = B(a,B( / 9, #(*))) (x), V$(f,<)e?[. (53) 
Differentiating (|52j) with respect to the parameter a, we obtain for a = 



r/ 



r/ 



C(tf («))(£) = — B(a, *(*))(£) =— U„(t,exp{ab}UZ l {t, tt(t))) (x) . (54) 
da a=o da s J ' a=o 



The operator C(-) determined by ()54|) is a generator of the one-parametric family of symmetry 
operators (|52|) . 

Note that the operator C(-) is not a symmetry operator for equation (|27|) since the parame- 
ters £ in the evolution operator U^t, •) (|42j) depend on a. Indeed, the parameters (<£) found 
from equation (|43|) 

fl(*» £ )| t=0 = (exp{ab}0|0|exp{ab}</»), 0(5) = l^ 1 (i, tf(t)) (f ), ^{x,t) £ 

include the parameter a explicitly. Therefore, equation (|54|) includes the derivatives of the 
evolution operator U^t, •) with respect to the parameters £, and ()54j) is different in form from 
the symmetry operator 1)51(1 . 

Nonlinear superposition 

The nonlinear superposition principle for the GPE ()27j) can be formulated in terms of the 
evolution operator 
Let 

®i(x,t) = Ux(t,ih)@), ^2{x,t) = U„(t,ip 2 )(x) G V\ (55) 

be two partial solutions to the GPE ()27J) corresponding to the initial functions ipi (x) , ip 2 (%) 
£ Vft , respectively. 
Then, the function 

V(X, t) = U„(t, dtp! + C 2 1p 2 ) (Xj 

is a solution of equation 1)27)1 which corresponds to the initial function ciipi(x) + c 2 ip 2 (x), 
ci,C2 G K 1 . Therefore, 

^{x,t) = U^t 1 c l U~\t 1 ^ l {t)) + c 2 liz\t^ 2 {t))){x) (56) 

is a solution of equation 1)27(1 which corresponds to the solutions ^>i(x, t), ^f 2 (x,t) of the 
form and equation ((o^|) is the superposition principle for the GPE (J27J). 

5 Examples 

3D case 

Consider equations @ and (j3J) in a 3-dimensional space with operators H.(t), V(t, \I/) of the form 
1 e .A 2 ,^t,.s - fc^ 



n(t) = — [p--A(x,t)\ -e(E(t),x) + -x z , (57) 
2m V c / 2 

t&, t) = V(x - y) = V (l - - yf) . (58) 
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The external field in the linear operator (|57|) is the superposition of a constant magnetic field 
H = (0,0, H) with vector potential A = x x, an electric field E(t) = (E cos Lot, E sin u>t, 0) 
periodic in time with frequency to, and the field of an isotropic oscillator with potential \x 2 , 
A; > 0. The operator H.(t) is the same as in |15j . and V{x — y) is obtained from the Taylor 
expansion of the Gaussian potential V(x — y) = Voexp — 2 $ • ^ n notations (|27 |) -(|3U |) . we 
have 

ftzz = ( V? p V, px ) , Ti-pp = — hx3, Hxx = diag{mw 2 ,mcj 2 ,mt^}; 

n z = (U p ,H x y, U p = 0, H x = -eE(t); (59) 

the nonzero elements of the 3x3 matrix 7i px (= T~i xp ) are: 7i. PlX2 = — W P2Xl = \loh- Here, 
u>H = ^ is a cyclotron frequency, lo ti i{v) 2 = is a "nonlinear frequency", lOq = ~, J\ = 
w o + (T)' " aVo'Mr,) 2 , J 2 = co 2 - CiVoMv) 2 , C(o) = sign(xa), a G R 1 , 77 = fy. The 
nonzero matrices in © are: W xa; = W yy = —W xy = —^VqI^xs. 

For the block matrix A(t) of the form (|40j) . the matrices Ai, A2, A3, A4 are of order 3x3 and 
they can be written as 

Ai = A 4 = -^-u(t), A 2 = -m 2 u(t), A 3 = -— r(t)u(i), 

where 

r . , N . , s . , ^ / cos(iu^) -sin^u^t) 

f (t = diag <^ v ; , v 1 ' , v M , u t = sin cos (±u H t) 

I wi wi cj 2 J V 1 

Here u(i) is an 50(3) matrix, uT(i)u(i) = u(i)u T (t) = I 3x3 . Denote by (P(t,h,£), X(t,h,€)) 
the general solution of the system (|16j) with (|59|). for which we use for short the notation 
(P(t),X(t)). We use similar notation for (jUJ, S(t,h,g(t,<t)) = S(t). 
Let us introduce the functions 



0(Ax J Ai/,t J a l P(t),P(a),w,W2 f) C) = J - - . -- exp I -(P(t)Ax - P{s)Ay) 

y 2nrnsm{co(t — sjj 

x ex p( o fe ■ \\ - s))(Arr 2 + Ay 2 ) - 2cos(^(t - a)) Ax Ay) 1 (60) 

1 2nsm(u]{t — s)) V 2 /) 

and 

G(Ax, Ay, t, s, P(t),P(s),co, C) = 0(Ax, Ay, t, s, P(t), P(s), to, u H , £)U=o- (61) 

Then, the Green function (|3*5|) reads 

y, t, s, g(t, £), g(s, £)) = £(Aa;i, Ayi, t, s, P\(t), Pi(s),ui,u> H , £) 
x C?(Aa; 2 ,Ay 2 ,t,s,P 2 (t),P 2 (s),u;i,u; /f ,G:)g(Ax3,Ay3,t,s,P3(t),P3(s),c 1 J 2 ,C:) 
* row o/_m! f * / sin(^(t-s)) 



x exp |-[5(t) - 5(a)] j exp j- I -ma^ ^Vj^ _ g)) (AxiAy 2 - Ax 2 A yi )J | . (62) 

The nonlinear evolution operator (|42j) is constructed with the use of (|62|) in which £ are 
determined by (|43[). 
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ID case 

To demonstrate symmetry operators in explicit form consider the one-dimensional case of the 
GPE (23) following [IE]. 

The operators ()57|) and ()58|) for n = 1 take the form 

H(z,t) = 1 — eExcosut, V(z, w, t) = - [ax 2 + 2bxy + cy 2 ] . (63) 

—Jill _ 

Here, x = x G R 1 , y = y £ R 1 ; k > 0, m, e, E, a, 6, and c are parameters of the potential. 

In notations we have fl(t, ft, g(t, £)) = ^P 2 (t, C)+ ^kX 2 (t, £)- e£X(i, £) cos wt+ 

f ca ra (t, £)+ f (a + 26 + c)X 2 (t, £), 



The Green function of the linear associated Schrodinger equation (|31|) follows directly from (|60(1 
and (EH). It reads 



GJx,y,t,s,g(t,<t),e(sX)) = G(Ax,Ay,t,s,P(t),P(s),n,<t}exp{ Us(t) - S(s)}} . (64) 



Here, uq = y/k/m and f2 2 = WQ+C(a)u; 2 ; (a). The evolution operator U^t, s, •) (|l2*|) is determined 
by (|64jl. where the constants £ are changed by determined by equation (|2*3*|) : g(s, <£) = 

0(C) = g^|g|^), 0(C) = (P(£),X(£), o-^C),^^), The LASE (El is 

known to have a partial solution in Gaussian form. For the operators (|63|) . this solution can be 
obtained as 



*<°> (x, t, fl (t, C)) = rf^e-*** exp{ | (5(t, ft, fl(t, C))+P(t, C) Ax) - ^ Ax 2 } . 

Setting tP(x) = $ 0) (x,0,g(0, <£)) = <4 0) (x,0) in & we can find the parameters <£ in the 
form C(y>) = Co = £(^ 0) (0)) = (V,x ,0,0, ^Y, where O 2 = u 2 + C(a + b)u 2 nl {a + 6), 

po = ^fm^/W^ff, an d x = (^ 0) |x|d. 0) )/ll^ 0) || 2 . Then, 

*$\x,t,£ ) = 1> 0) (x,t,£i(t,C(V))) (65) 

is a partial solution of the 1-dimensional GPE H27|) . 

Let us change the operators a in relation (|51|) by the operators a + and a of the form 

: [Aj}q — imf2Ax ] , a + = [Apo + imQAxo] , 



where Ap = —ihd x — po and Ax = x — xq. Then, the operators A^(-) determined by the 
relations 

(*(*)) (x) = a+ f/" 1 (t, *(*))) (x), (*(t)) (x) = a C/- 1 (i, *(*))) (x), 

where ^(x,i) E "P^, are the symmetry operators of equation (|27|). With the operators 1(±)(-), 
we obtain the set of solutions for the GPE 



*(°)(x,t, (t,£ n )) = -±={A(+\.)) n ^\x,t,Q(t,£ )) 



Exact Solutions and Symmetry Operators 
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inUt( 1 Yrr fJ mn A„\,TM 



where n G Z + , H n (£) are Hermite polynomials and 
-T/,s Z' eE h(2n+l)\ J 

It can be verified that 

1« (* (?) (t, fl (t, e n ))) (x) = v^TT ^(s, t, fl (t, e n+1 )), 
l(-)(*W(t,g(t,e:n)))(x) = v^<°2i(^t,0(t,^i)) 

and the operators A^(-) are nonlinear analogs of "creation-annihilation" operators. It can be 
noted that the functions ^n\x, t, g(t, satisfy the quasi-periodic condition 



¥ 0) (x t + T £ T ) = e - i£nT/h ^ 0) (x t £ q 
where £ n is the quasi-energy. The spectrum of quasi-energies is given by the relation 
e 2 E 2 e 2 E 2 [oj 2 -uj 2 - C(a + 2b + c)u^(a + 2b + c)] 



2m(Q 2 -LU 2 ) 4m(Q 2 -uj 2 ) 2 

)(' 



+ h(tt + -^tt) (n+ - 
V 2mQ,)\ 2 

6 Concluding remarks 

The evolution operator (|42[) obtained in Sectional in explicit form enables one to look in a new 
fashion at the problem of construction of semiclassically concentrated solutions to the nonlocal 
Gross-Pitaevskii equation (J2J). In particular, the semiclassical asymptotics constructed in jl3M 14j 
can be presented in more compact and visual form. In addition, the evolution operator leads to 
the nonlinear superposition principle and to the general form of symmetry operators (|51[) . 
The latter can be used for study of symmetries of the Gross-Pitaevskii equation under con- 
sideration as long as direct finding of the symmetries by means of solution of the determining 
equations is a nontrivial problem. 
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